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ABSTRACT 



Aims. The observed brightness of Type la supernovae is affected by gravitational lensing caused by the mass distribu- 
tion along the line of sight, which introduces an additional dispersion into the Hubble diagram. We look for evidence 
of lensing in the SuperNova Legacy Survey 3-year data set. 

Methods. We investigate the correlation between the residuals from the Hubble diagram and the gravitational magni- 
fication based on a modeling of the mass distribution of foreground galaxies. A deep photometric catalog, photometric 
redshifts, and well established mass luminosity relations are used. 

Results. We find evidence of a lensing signal with a 2.3 a significance. The current result is limited by the number of 
SNe, their redshift distribution, and the other sources of scatter in the Hubble diagram. Separating the galaxy pop- 
ulation into a red and a blue sample has a positive impact on the significance of the signal detection. On the other 
hand, increasing the depth of the galaxy catalog, the precision of photometric redshifts or reducing the scatter in the 
mass luminosity relations have little effect. We show that for the full SuperNova Legacy Survey sample (~400 spectro- 
scopically confirmed Type la SNe and ~200 photometrically identified Type la SNe), there is an 80% probability of 
detecting the lensing signal with a 3 a significance. 

Key words, supernovae: general - cosmology: observations - gravitational lensing 



* Based on observations obtained with 

MegaPrime/MegaCam, a joint project of CFHT and 
CEA/DAPNIA, at the Canada- France-Hawaii Telescope 
(CFHT) which is operated by the National Research Council 
(NRC) of Canada, the Institut National des Sciences de 
l'Univers of the Centre National de la Recherche Scientifique 
(CNRS) of France, and the University of Hawaii. This work 
is based in part on data products produced at the Canadian 
Astronomy Data Centre as part of the Canada-France-Hawaii 
Telescope Legacy Survey, a collaborative project of NRC 
and CNRS. Based on observations obtained at the European 
Southern Observatory using the Very Large Telescope on the 
Cerro Paranal (ESO Large Program 171.A-0486 & 176. A- 
0589). Based on observations (programs GS-2003B-Q-8, 
GN-2003B-Q-9, GS-2004A-Q-11, GN-2004A-Q-19, GS-2004B- 
Q-31, GN-2004B-Q-16, GS-2005A-Q-11, GN-2005A-Q-11, 
GS-2005B-Q-6, GN-2005B-Q-7, GN-2006A-Q-7, GN-2006B-Q- 
10, GN-2007A-Q-8) obtained at the Gemini Observatory, which 



1. Introduction 

Type la supernovae (SNe la) have become an essential tool 
of observational cosmology (Riess et al. 1998 Perlmutter 
[eTaL] p99l |Astier et ed] |2006| |Wood-Vasey et al.| |2007| 



is operated by the Association of Universities for Research 
in Astronomy, Inc., under a cooperative agreement with 
the NSF on behalf of the Gemini partnership: the National 
Science Foundation (United States), the Particle Physics 
and Astronomy Research Council (United Kingdom), the 
National Research Council (Canada), CONICYT (Chile), the 
Australian Research Council (Australia), CNPq (Brazil) and 
CONICET (Argentina). Based on observations obtained at 
the W.M. Keck Observatory, which is operated as a scientific 
partnership among the California Institute of Technology, the 
University of California and the National Aeronautics and 
Space Administration. The Observatory was made possible by 
the generous financial support of the W.M. Keck Foundation. 
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Riess et al.j [20071 |Freedman eTaT] [20091 |Kessler et all 
2009 ) . By studying the distance-redshift relation of a large 



number of SNe la over a wide range in redshift, the equation 
of state of dark energy can be measured. Although SNe la 
can be calibrated to be good standard candles, they are 
affected by systematic effects such as extinction by dust 
and gravitational lensing. Because of gravitational lensing, 
most supernovae will be slightly demagnified and some will 
be significantly magnified due to the mass inhomogeneities 
along the line of sight. 

This causes an additional dispersion in the observed 
magnitudes and thus in the Hubble diagram. The effect 
is greater at high redshift, but is still less than the scatter 
in the Hubble diagram for a redshift up to 1, as we see in 
this paper. 

Magnification of SNe la can be estimated in 2 ways. The 
Hubble diagram residuals, computed assuming the best fit 
cosmological model, give 

an indirect measure of the SN magnification. The mag- 
nification can, on the other hand, also be estimated by 

modeling the foreground galaxy mass distribution us- 
ing galaxy photometric measurements together with de- 
rived mass-luminosity relations for galaxies and dark mat- 
ter halo models, ff a correlation between these two esti- 
mates is found, it could give interesting insight in the mod- 
eling of the mass distribution of the foreg round galaxies 



This will be discussed in another paper ( |Jonsson et al. 
20091 



The lensing si gnature in supernovae sam ples has already 
been sought for 



Williams & Song 
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with chosen mass-luminosity relations) and the residuals 
from the Hubble diagram. 

Evaluating the expected magnification by foreground 
galaxies can be schematically split into three steps : first ob- 
taining redshifts and absolute magnitudes of the foreground 
galaxies, then converting these informations into mass, and 
finally turning this set of mass and redshift estimates into 
a magnification estimate. 

The outline of the paper is as follows, ^introduces the 
mass-luminosity relations used in this analysis. In f|3 the 
SNLS 3-year data set is presented, ^describes the ana ysis 
steps to estimate the supernovae magnifications, while in f|5] 
we present the measured correlation and compare it with 
expectations from simulations. Finally, we summarize in 

2. Mass Luminosity relations 

There are several ways of inferring the combined mass of a 
galaxy and its dark matter halo using the measured lumi- 
nosity. We make use of two approaches in this paper. 

i) The line-of-sight velocity dispersions in elliptical 
galaxies and rotational velocities in spiral galaxies are corre- 
lated with their luminosities. Those are the well established 
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Faber & Jackson 
77 ' Haynes et al. 



Search Team and the Supernova Cosmology Project with 
the density of the foreground galaxies. They found that 
brighter supernovae preferentially lay behind over-dense re- 
gions at a 99 % CL. However, they find a lensing contribu- 
tion to the scatter in the Hubble diagram of 0.11 mag (for 
a median redshift below 0.5) which is larger than expected 
(for example fro m lensing statistic s ), as noted by the au- 
thors themselves. |Menard &: Dalai ( |2005[ ) remark that the 
measured size of the effect is also incompatible with shear 
variance measurements. Finally, such a large magnification 
effect would cause a significant increase o f Hubb l e diag ram 
scatter with redshift, undetected today. Wang (2005) de- 
rived the expected weak lensing signatures of Type la SNe 
by convolving the intrinsic distribution in peak luminos- 
ity with expected magnification distributions and compared 
the expected and measured skewnesses of the residual dis- 
tribution of 110 hi gh and low-z SNe from the Riess sample 
( Riess et al. 2004[ ) . The signal found is dominated by high 
redshift events uncorrected for the NICMOS non-linearity. 
Its origin remains unknown and its significance is not pro- 
vided. Partl y using the same supernova sample, Menard & 
Dalai (2005 1 searched for a correlation using SDSS photom- 
etry of the foreground galax ies and found no evid ence for 
a correlation. More recently Jonsson et al. (2007) found a 
weak correlation with a confidence level of 90% using high- 
z supernovae from the GOODS fields. A firm detection of 
this effect therefore remains to be found. 

The SNLS sample is currently the best suited sample 
for a possible detection of such a si gnal thanks to its lar ge 
number of SNe la at high redshift ( |Jonsson et al.j 2008). 

In this paper, we investigate the possible correlation be- 
tween the magnification of the SNLS SNe derived from pho- 
tometric measurements of the foreground galaxies (together 



Faber- Jackson (Poveda 
[1976] and Tully-Fisher 
1999) relations . 

ii) Galaxy-galaxy lensing measurements allow a relation 
to be determined between luminosity and mass once a func- 
tional form is assumed for the mass density profile. 

We present below recent results for these approaches 
and compare them assuming a Singular Isothermal Sphere 
dark matter density profile (hereaft er SIS) which has prove n 



to be a good fit to lensing galaxies (Koopmans et al. 



2006) 



2.1. Tully-Fisher relation for spiral galaxies 



We have chosen to use the Boehm et al. (2004) results 
for the Tully-Fisher (TF) relation. The maximum rota- 
tional velocities of 77 spiral galaxies in a redshift range 
of 0.1 < z < 1.0 were measured in the FORS Deep Field. 
Anchoring the TF rela tion at low z using the results from 
Pierce & Tully (19921, they obtain the following relation 
between the maximum rotation velocity w max and the abso- 
lute magnitude Mb of the galaxy in the rest-frame B-band 
(in a Vega magnitude system): 

log w max = -0.134 (M B + (3.61 ± 0.24) + (1.22 ± 0.56) • z) 

(1) 

The observed scatter of Mb about this relation is 0.41 mag 
(r.m.s.). The redshift dependence expresses the fact that a 



a younger stellar population (|Barden et al.l 2003; Milvang- 


Jensen et al. 

1 TTTT - — ' 


2003 


Bamford et al. 2006 Chiu et al. 


2008). 



ies, measured at large galactic radii, is dominated by the 
dark matter mass. Assuming a SIS dark matter density 
profile for which the rotational velocity is constant, we 
have w max — Vmt{r — > oo) = \/2er, where o is the one- 
dimensional r.m.s. of the velocity distribution. 

2.2. Faber- Jackson relation for elliptical galaxies 



We ad opt the Faber-Jackson relation (FJ) by Mitchell et al. 
(|2005[) (hereafter M05) derived for a sample of ~ 30, 000 el- 
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liptical galaxies using SDSS data. The selection criteria for 
the sample and the estimate of the velocity dispersion are 
described in Bernardi et al. (20031. The observed velocity 
dispersion has been determined by analyzing the integrated 
spectrum of the whole galaxy and aperture corrected to a 
standard effective radius. M05 find the following relation 
between the velocity dispersion and the absolute magni- 
tude M r of the galaxy in the rest-frame r-band (in a AB 
system) : 

< log(er) >= 2.2 - 0.091(M r + 20.79 + 0.85z) (2) 

As for the TF relation, the redshift dependence accounts 
for the evolution of the age and hence luminosity of the 
stellar population. The scatter in the FJ relation induces 
an uncertainty in the e stimate of the velo city dispersion 
which has been given by Sheth et al. ( 2003 ) : 



r.m.s.(logcr) = 0.79(1 + 0.17(M r + 21.025 + 0.85z)) 

The measured velocity dispersions are aperture-corrected 
central velocity dispersions which are almost equal to the 



dark matter velocity dispersions (see Franx 1993 



Kochanek 



19941, so we can identify those measurements with the ve 



locity dispersion parameter of the SIS profile. 

We choose to translate the FJ r-band relation into a 
restframe B-band relation since restframe B is covered with 
SNLS griz bands up to z ~ 1. To convert the SDSS r-band 
absolute magnitudes in the AB system to standard B-band 
Vega absolute magnitudes, a typical color Mb — M r for 
ellipticals in the AB system is estim ated yielding M B — 
M r = 1.2CQ ( |Gunnarsson et al.| |2006[ ) and an AB to Vega 
relation Bab = -#Ve ga — 0.12 is adopted. 



where L is the luminosity of the galaxy. The fiducial lu- 
minosity, L* = 10 10 /i -2 L r o is given in the SDSS r-band. 
For a conversion to the -B-band, they have calculated that 
the galaxies in their sample with a fiducial luminosity of 
L* = 10 10 Ii~ 2 Lbq in the B-band have a fiducial luminos- 
ity of L» = 1.1 x W 10 h- 2 L r(3 in the SDSS r-band. Note 
that for the SIS profile, the halo is characterized by its ve- 
locity dispersion which is the key parameter ente ring in the 
estimation of the magnification (see section 4.4 ) . 



2.4. Comparison assuming a Singular Isothermal Sphere 
profile 

In Figure [T] we show the velocity dispersion as a function 
of absolute B-band magni tude fo r the three different mass- 
luminosity relations (Eq. |1|2|3 for TF, FJ relations and 
galaxy-galaxy lensing results of K06) assuming a SIS dark 
matter mass profile. The mean redshift of the lens galaxies 
in the K06-relation is ~ 0.4. For a given B-band luminosity, 
elliptical galaxies are more massive than spirals so the FJ 
relation naturally gives higher velocity dispersion than the 
TF relation. The K06 relation and the FJ relation appear 
quite similar. For high luminosity galaxies the difference in 
mass estimate from the three relations can lead to very dif- 
ferent magnifications. As an example, for a bright galaxy of 
absolute magnitude Mg—-22, the velocity dispersions range 
from 170 km s" 1 (TF) to 249 km s" 1 (FJ) and 248 km s" 1 
(K06). As a consequence, in order to obtain accurate es- 
timates of the magnification it is important to take into 
account the galaxy type. 



2.3. Galaxy-galaxy lensing mass estimates 

The galaxy-galaxy lensing signal manifests itself by im- 
ages of background (faint) g alaxies being distorted by fore 
ground (brighter) galaxies (Hoekstra et al.[ [20041 [2005 



Parker et ahj |2005[ |Kleinheinrich et al.[|2006|lMandelbauml 
et al.[ 120061). Unfortunately, one can only study ensemble 



averaged properties because the weak-lensing distortion in- 
duced by an individual galaxy is too small to be detected. 
The measured quantity is the mean tangential shear which 
is proportional to the total lens galaxy-halo mass. Since the 
lensing signal depends on the angular diameter distance be- 
tween observer, lens and source, galaxy-galaxy lensing re- 
sults strongly depend on the properties of the galaxy sam- 
ple, which differ for different surveys. 

We use the ga laxy-g alaxy lensing result from 



Kleinheinrich et all ( |2006| (hereafter K06). The data 
is taken from the COMBO-17 survey which consists of 
observations in five broad-band filters (UBVRI) and 12 
medium-band filters. They report a mass to luminosity 
relation for the full sample probed out to a maximum 



radius of 150 
profiles: 



h 1 kpc when modeling the lenses as SIS 



(T 
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+ 18 
24 



L 



10 10 h~ 2 L 



r& 



n ,„+0.12 



km.s 



1 With the galaxy spectral sequence described in [4.2 
for the range of rest-frame U 



(3) 



and 



V colors we consider to select 
photometrically elliptical galaxies, we obtain on average Mb — 
M r — 0.86 with an r.m.s of 0.07. This corresponds to an average 
change of a of 7% with an r.m.s of 1% (see Eq. [2j|. We ignore 
this small correction in the following. 




Fig. 1. Velocity dispersions as a function of the abso- 
lute magnitude in the B-band for three different mass- 
luminosity relations at a redshift of z = 0.4. In green (dash- 
dotted) the K06 relation, in red (dotted) the FJ relation and 
in black (solid) the TF relation. The plot shows the velocity 
dispersion in km s _1 as a function of absolute magnitude 
in the B-band. 
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3. Data 



4.1. Galaxy photometric catalog 



The SuperNova Legacy Survey consists of an imaging sur- 
vey which is part of the (deep component of the) Canada- 
France-Hawaii Telescope Legacy Survejj^] and a spectro- 
scopic survey done on 8m class telescopes. The imaging was 
done with the MegaCam imager (360 Megapixels, 1 deg 2 ) 
which detects and monitors the light curves of the SNe, 
over four different fields of 1 deg 2 (called DL D2, D3, D4). 
Dl overlaps partly with the VVDS survejj^J D2 with the 
COSMOS survejQ and D3 with the DEEP2 survejQ We 
will later use the spectroscopic galaxy redshifts from the 
DEEP2 and VVDS surveys to establish and test photomet- 
ric galaxy redshifts. 

The "rolling" search method was used for the griz 
bands, which means observing the same field every 3-4 days 
during dark and gray time (typically for a period of ~ 16 
days around new moon) for as long as it remains visible 
(from 5 to 7 months a year). Images in u are also used 
in this analysis. The images are used both to monitor the 
SN light curves and to build photometric catalogs of the 
objects in the fields, in particular the galaxies. 

For each field, the total exposure times used to build 
the galaxy catalogs for this analysis amount to more than 
20 h per band, reaching typically 60 hours in the i-band. 

Spectroscopy of the supernovae is crucial in order to 
obtain redshifts, and to confirm the type of each super- 
nova. The spectroscopic follow-up for the SNLS was done 



with the VLT (Baumont et al. 2008 Balland et al. 2009) 



Gemini ( Howell et al. 
telescopes (Ellis et al 



. i 200 5]]Bronder et al.[|2008j and~geck 
L , I2008D- The 3-year data set used in 



this paper consists of 233 spectroscopically confirmed Type 
la supernovae used for cosmological analyses. A detailed de- 
scription of the S N survey and the SN data ana l ysis m eth- 
ods is provided in |Astier et"aLl ( |2006[ ); |Guy et al.| ( 2009[ ). An 
important feature, relevant here, is that the SNLS Hubble 
diagram exhibits a r.m.s scatter of 0.16 mag. 



4. Estimating the supernovae magnification 

We describe in this section the analysis steps followed to 
calculate the supernovae magnifications. We need a cata- 
log of the field galaxies with an estimate of the redshift 
and the rest-frame B, V and U band absolute magnitudes 
for each galaxy. Computing the restframe magnitudes re- 
quires, in addition to the redshift, the knowledge of the 
galaxy SED. The rest-frame £?-band absolute magnitude is 
used to convert the luminosity of the galaxy into a mass 
estimate using the mass-luminosity relations introduced in 
section [2] whereas the U—V color is used to separate the 
galaxies into a red and a blue population. The photometric 
catalog is described in [4.1 The method we have used to 
determine photometric redshifts and absolute magnitudes 
is accounted for in |4.2| t he galaxy photometric classifi- 



cation is presented in [4.3 and the estimation of the SN 



magnification is done in §4.4| 



2 see http://www.cfht.hawaii.edu/Science/CFHLS 

3 http:/ /www. oamp.fr/virmos/vvds. htm 

4 http:/ /cosmos. astro. caltech.edu/ 

5 http://deep.berkeley.edu/ 



The galaxy catalogs are built on deep image stacks in the 
ugriz filters. The deep stacks are constructed by selecting 
80 % of the best quality images (6241 images for the four 
fields). Transmission and seeing cuts (e.g. fwhm < 1.15") 
are applied. Because we have fewer exposures in the it-band 
than in the others, less stringent quality cuts are applied 
to these images. The selected images are co-added using 
the SWARP v2.10 package]^] The source detection and pho- 
tometry i s perfo rmed using SExtractor V2.4.4 ((Bertin & 
19961) ) in double image mode. The detection is 
i ban 



Arnouts 
made in 



the i band. 



A cut on the signal-to-noise ratio in the i-band, S/N > 
15 has been made so as to maximize the signal detec- 
tion while minimizing the many spurious detections around 
stars halos. We then use the AUTO SExtractor flux, com- 
puted in an elliptic aperture to extract the galaxies. The 
different cuts lead to a limiting magnitude i = 25 (Vega 
magnitude system). Zero points are computed using cir- 
cular photometry on a tertiary stars catalog, described in 
(Regnault et al. 



2009) 



Two categories of objects have to be identified in the 
catalog : stars, and the host galaxies of the SNe. In order 
to identify stars, we estimate the second moments of all 
objects in the catalog using a 2D Gaussian fit. We then 
look for an accumulation in the second moment m xx — m yy 
space, taking into account the variations of the PSF across 
the focal plane. Objects in this clump are assumed to be 
stars. 

Since the photo-z of the host galaxy is less precise than 
the spectroscopic z of the SN it is possible to have the con- 
figuration of the host galaxy being in the line-of-sight of 
the supernova and hence wrongly contributing to its mag- 
nification. As a consequence it is important to identify and 
exclude the host galaxy from the analysis. Note that identi- 
fying the host galaxy on the images is not always obvious. 
The host galaxy is identified using two criteria: the min- 
imum distance to the supernova location and the match 
between the photometric redshift of the galaxy and the 
spectroscopic redshift of the supernova. 

To estimate these quantities it is necessary to use im- 
ages uncontaminated by the supernova light. For a given 
supernova, we exclude the images taken during the same 
season, i.e. the 6 consecutive months during which the field 
was observed. The photometric redshift is derive d u sing 
the SNLS photometric redshift code (see section 4.2 1. A 
normalized distance d is computed using SExtractor shape 
parameters, so that d < 1 defines the photometric ellipse of 
the galaxy. When no object is found within d < 1.3 of the 
SN, the host of the SN is undetected. When more than one 
object is detected close to the supernova location, we check 
the match between the galaxies photometric and the super- 
nova spectroscopic redshifts. Ambiguous cases are nagged 
as problematic and can lead to exclusion of the SN if the 
uncertainty in the determination of the host galaxy has an 
important impact on the magnification of the SN in ques- 
tion. Four SNe from the sample were initially excluded in 
this way. 

However, we are fortunate to have HST imaging from 
the COSMOS field (IScoville et all [20071 IKoekemoer et al ' 



2007) together with newly published high resolution red- 



http://terapix.iap.fr/soft/swarp/ 
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shifts ( jllbert et al. |2009[ ), also from the COSMOS for galax- 
ies in the D2 held. As a consequence, two of the initially 
excluded SNe have been reexamined, SNLS-04D2kr and 
SNLS-05D2bt. These two SNe are detected very close to 
a presumed host galaxy which photometric redshift does 
not match the spectroscopic redshift of the SN. Figure [2] 
and E] show the SN location on the CFHT and the HST 
images for the two SNe in question. 

For SNLS-04D2kr at z = 0.744, a smaller and hardly 
visible galaxy is detected at the location of the SN in the 
HST image, very close to the presumed - large - host 
galaxy. The r edshif t assigned to this la rge galaxy from 



( Ilbert et al] [2006 1, ( |Ilbert et al^ |2009| ), and the SNLS 



photometric redshift code is z = 0.168, 0.228 and 0.3 re- 
spectively, implying that the large galaxy in question is not 
the host galaxy, but a foreground galaxy. The host galaxy 
is likely to be the small elongated galaxy detected in the 
HST image. 

The same situation arises for SNLS-05D2bt at 
z = 0.68 : the presumed host galaxy has a photometric 
redshift of z = 0. 31 and 0.32 from t he SNLS photometric 
redshift code and ( Ilbert et al. 2006 1 respectively. 



We obtain in this way a photometric catalog in the ugriz 
bands, where stars and supernovae host galaxies are iden- 
tified. The host galaxy magnitudes are free from supernova 
light. 



4.1.1. Masks 

Halos and diffraction spikes around bright stars generate 
spurious galaxy detections in the catalog. The size of the 
halos and the shape of the spikes are fixed by the geome- 
try of the telescope optics, but their intensity in the images 
is directly proportional to the brightness of the stars. Also, 
for the brightest stars, pixels reaching their saturation level 
induce bleedings in the CCD images which have to be ac- 
counted for. 

Due to those effects, circles with radius varying from 
50 to 600 pixels (10 to 120") centered around bright stars 
are masked out. Those masked regions represent 22% of the 
field of view. An example is displayed in Figure [4] for the 
Dl field. 

For each superno va w e have included galaxies within 
a radius of 60" (see £4.4). This will be referred to as the 



selected area. In theory, one should exclude all SNe where 
the selected area overlap a masked area, but this means 
rejecting about half of the SNe. However, the effect on the 
magnification of small spurious galaxies or diffraction spikes 
in the outskirts of the selected area is small and we chose to 
keep all SNe with an angular separation larger than 40"of 
masked areas. This leads to the exclusion of 53 SNe. 

All SNe have also been looked at by eye to see whether 
there are other effects that could lead to exclusion, e.g. 
non-masked stars very close to the line of sight leading to 
the possibility of excluding an important galaxy hidden be- 
hind the star or leading to possible bad photometry for the 
surrounding galaxies. Looking at the SNe we exclude 9 of 
them and as a conclusion we keep 171 SNe out of 233 in 
the initial sample. 



♦ • 
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Fig. 4. A section of the Dl field (5000x4500 pixels or 
15.4x13.9') with the masked areas shown as purple circles. 



4.2. Photometric redshifts 
High quality phot ometric redshifts have been published by 



Ilbert et al. ( 2006 ) for the galaxies in the SNLS fields down 
to zab — 24. In these catalogs, fluxes of galaxies hosting su- 
pernovae are affected by the supernovae light, leading some- 
times to unreliable typing and redshift for these galaxies. 
We also need the SED used to derive the photometric red- 
shift, so as to compute the restframe magnitudes. Finally, 
we need to be able to propagate easily the uncertainties of 
the photometric measurements to the photo-z and the ab- 
solute magnitude estimations. For these reasons, and so as 
to control the error propagation path, we have chosen to de- 
rive the photometric redshifts and the absolute magnitudes 
as follows. 

We first define a continuous one parameter (a*) galaxy 
spectral sequence, F(a*,A). Indeed as we have two other 
parameters to estimate (magnitude and redshift), we can- 
not afford more than a single parameter to index the diver- 
sity of galaxies in order to keep one constraint, when, even 
though 5 bands are used, only 4 are well measured {griz, u 
being significantly shallower) N To define this spectral tem- 
plate, we use the galaxy evo lution model PEGASE.2 (|Fioc| 



1999[ ) , using a variety of galaxy 



& Rocca-Volmerange 

law oc (t/r) x exp(— t/r) and galaxy ages. The initially zero 
gas metallicity evolves through the successive generations 
of stars and ranges from Z~0. to 0.03. Extinction is com- 
puted using a transfer model for an inclination-averaged 
disk distribution, the optical depth being estimated from 
the mass of gas and the metallicity. We left aside the pos- 
sibility to add an additional extinction. 

The next step consists in optimizing the spectral se- 
quence so as to reproduce the observed colors of our data 
in the best way. 

It can be corrected to about 30 %, so that the final 
template spectra do differ from the initial computed tem- 
plate : our technique is thus marginally sensitive to the in- 
put model. The training set comprises a sample of galaxies 



7 It would be conceivable to add an extinction parameter if 
we restricted our scope to bright galaxies. 
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(a) CFHT image. 



(a) CFHT image 





(b) HST image from the COSMOS survey 



(b) HST image from the COSMOS survey 



Fig. 2. Images of the field around the supernova SNLS- 
04D2kr at z = 0.744. The images span 30" and are ori- 
ented north- up/east-left. A red circle indicates the super- 
nova position. A blue ellipse indicates the closest galaxy, 
as detected on the CFHT image. 



Fig. 3. Field around SNLS-05D2bt at z = 0.68. As for 
SNLS-04D2kr, a smaller and fainter galaxy is conspicuous 
on the SN location on the HST image. 



with known spectroscopi c redshift from the DEEP2 survey 
( Davis et al. , 2003 , 20071) . The spectral sequence consists in 
the galaxy flux as a function of wavelength and the age in- 
dex, and we optimize a multiplicative correction (a smooth 
function of wavelength and age index), along with offsets 
to the photometric zero points (which can be interpreted as 
differential aperture corrections). We find a residual scat- 
ter of 0.097, 0.015, 0.035, 0.025, 0.050 mag for ugriz bands 
respectively. 



The performance of the photometric redshift compu- 
tation is evaluated using VVDS spectroscopic redshifts 
avail able for the Dl field (3 595 galaxies at 0.01 < z < 



1.5) (Le Fevre et al. 2004 1. The redshift residuals (i.e. 



Az=photometric redshift - spectroscopic redshift) as a 
function of spectroscopic redshift are shown in Figure [5j 
For «ab < 24, the fraction of catastrophic failures for which 
Az/(1 + z) > 0.15 is 6%. Eliminating catastrophic failures, 
we obtain <j& z — 0.065 and cta 



z/(l+z) 



= 0.037. Note that 



we do not apply any prior on the photometric redshift, 



T. Kronborg et al., SNLS collaboration: Gravitational lensing of the SNLS sample. 



7 



as a large fraction of the SNLS galaxies are much fainter 
than those of the training set. 
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0.5 1 
spectroscopic z 

Fig. 5. Performance of the photometric redshift determi- 
nation: the difference Az between spectroscopic and pho- 
tometric redshift as a function of spectroscopic redshift. 
The sp ectroscopic redshifts w ere obtained from the VVDS 
survey (Le Fevre et al. 2004). 



Uncertainties are estimated using a Monte Carlo prop- 
agation of magnitude uncertainties accounting for the mea- 
surement uncertainties and the residual scatter obtained 
from the training procedure. In Figure [6] we compare the 
scatter of Az for the VVDS test sample to the uncertainty 
on the photometric redshift derived from the Monte Carlo 
propagation. Both are in reasonable agreement which vali- 
dates our method for propagating uncertainties. 
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Fig. 6. Performance of our photometric redshift determi- 
nation: the uncertainty in the estimated photo-z using the 
spectroscopic VVDS redshift as a function of photo-z (in 
dotted red) and the photo-z uncertainty (solid black line) 



For some of the galaxies it has been possible to ob- 
tain spectroscopic redshifts from VVDS, DEEP-2 and the 
SNLS spectroscopic program on VLT (some field galaxies 
were targeted at the same time as the SN primary target 



during observations with FORS-2 in multi-slit mode). The 
D2 field overlaps with the COSMOS field and it has thus 
been possible to use high resolution photometric re dshifts 
for a large fraction of the galaxies in the D2 field (Ilbert 
etlll|2009). 



4.3. Classification of spiral and elliptical galaxies based on 
their colors 

The Tully-Fisher and Faber- Jackson relations are derived 
for spiral galaxies and elliptical galaxies respectively and as 
a consequence it is necessary to separate the SNLS galaxies 
into spirals and ellipticals. Note that this classification is 
not strictly necessary: we could blindly apply some sort of 
average mass-luminosity relation to all galaxies. There is 
indeed no typing involved in our analysis based on galaxy- 
galaxy lensing mass estimates. Separating in broad types is 
not more than a way to a priori improve the significance of 
a potential detection. 

The SNLS imaging data does not permit a good mor- 
phological classification and we resorted to using a sim- 
ple color cut to classify our galaxies. The rest-frame color 
U — V is computed giving rise to two well separated distri- 
butions (the red and the blue population). In Figure [7] the 
rest-frame color U—V for the SNLS galaxies is shown. For 
U—V > 0.54, the galaxy is classified as an elliptical galaxy 
or else it is classified as a spiral galaxy. Using the rest-frame 
U—V color as a proxy for galaxy classification is obviously 
not optimal, and the effect of miss- classif ying the galaxies 
will degrade the lensing signal (see [ 5.3.3). 



4000 



2000 




U-V 

Fig. 7. Distribution of the rest-frame U—V color for the 
SNLS galaxies. A cut at U — V = 0.54 separates the distri- 
bution into red (elliptical) and blue (spiral) galaxies. 



To evaluate whether this cut provides a good estimate 
for galaxy t yping we compare o ur g alaxy types to thos e 
~~ p006| and lllbert et all p009} 



published in Ilbert et al. 



They correspond to the galaxy template type, as fitted with 
eventually some additional extinction, using respectively 5 
and 30 band magnitude measurements. For iab ^ 24, we 
have classified 89.5% of the COSMOS and 95.3% of the 
CFHTLS elliptical galaxies as ellipticals. However 56.8% 
and 26.9% of our elliptical galaxies are classified as spiral 
galaxies in the COSMOS and the CFHTLS respectively, 
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with an additional extinction of < E(B—V) >~ 0.2. In con- 
clusion, we identify mainly all elliptical galaxies, but there is 
a contamination of our color-selected elliptical sample with 
red extinct spiral galaxies. This contamination will dilute 
the lensing signal : the red spirals misidentified as ellipticals 
will be attributed too high a mass, falsely increasing the 
supernova expected gravitational magnification, and thus 
decreasing its correlation with Hubble diagram residuals. 

4.4. Gravitational magnification 

The density profile of a SIS can be written 



<7 

2ttt 2 



(4) 



The total mass of the halo, m(r) — 2a 2 r, diverges and as a 
consequence we can use a truncation radius, r t to be able to 
estimate the total mass of the halo. The truncation radius 
is chosen to be ^200 = V2<r/10 H{z) which is defined as the 
radius within which the mean mass density is 200 times the 
critical density. In this case, the convergence is given by 



ksis(0) 



7T0 



arctan < 



' 200 

6 2 D? 



(5) 



in the thin lens approximation, where 9 is the transverse 
angle of the image to the lens, and 6e = ^{c/c) 2 Di s /D 8 
is the Einstein radius. _D; S , D s and Di are respectively the 
angular distances from the lens to the source, from the ob- 
server to the source and from the observer to the lens, and 
cr is the velocity dispersion. In order to compute the mag- 
nification due to several ga laxy lenses, we use the publicly 
available software Q-LET (Gunnarsson 20041, which uses 
the multiple lens plane method. 



0.03 



^.02 



3 
b 



0.01 
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Fig. 8. r.m.s of the difference of magnifications n(d max ) — 
A*(<inax = 125") as a function of d max . This simulation (from 
real lines of sight) is done for a source at a redshift of 1. 



We have investigated the effect on the estimated mag- 
nification of truncating the galaxy list around lines of 
sight. We have considered for this purpose random lines 
of sight in the galaxy catalog, and computed for each of 
them the magnification for a SN at z = 1, considering 



sequentially an increasing number of galaxies in the fore- 
ground according to their (transverse) angle 9 to the SN. 
Figure [8] displays the r.m.s of the difference of magnifica- 
tions fi(d m „ < 125") — /i(^max = 125") as a function of 
dmax, where d max is the maximum transverse angle of galax- 
ies used to compute the magnification^ For d max = 60", the 
average error on the estimated magnification due to this 
galaxy selection is of 0.003. This number is much smaller 
than the other sources of uncertainties so that we can safely 
ignore galaxies at larger angles. 

4.4.1. Normalization of the magnification distribution 

To estimate the magnification using Q-LET, the lensing 
galaxies have been put on top of a homogeneously dis- 
tributed universe leading to computed gravitational magni- 
fications always greater than 1 compared to such a universe. 
As a consequence the magnifications should be normalized 
in order to force the average magnification to be equal to 
1, as imposed by flux conservation. 

For this purpose, we have considered numerous random 
lines of sight (randomly chosen source positions) using the 
true galaxy catalog for a range of source rcdshifts and com- 
puted the magnification for each of them. The average mag- 
nification for each redshift bin (of 0.1) is recorded and sub- 
sequently used to normalize the magnification estimates. 
This correction function is calculated for each field and each 
mass-luminosity relation. 

4.4.2. Magnification uncertainty calculation 

The uncertainties on the magnifications are evaluated with 
a Monte Carlo propagation of the magnitude uncertainties 
of the galaxies from the catalog and the scatter in the mass- 
luminosity relations. 

Magnitude uncertainties include both measurement un- 
certainties and the resi dual scatter of the photometric red- 
shift training 



sec 



4.2| 



A fit of the photometric redshift is performed for each 
Monte Carlo realization of each galaxy, except for those 
which have much more precise spectroscopic or photomet- 
ric redshifts from other sources. For instance, high resolu- 
tion photometric r edshifts are availa ble for galaxies of the 

For those, instead of 

fitting for the reds! 



ift. 



t et al.[|2009| ). 

we considered random realizations 



of the redshift within the errors given in their paper, and 
used our code with those fixed redshifts to determine the 
uncertainties on the absolute magnitudes. 

The scatter in the TF and FJ mass- lumi nosit y rela tion s 
are well defined from the observations (see ^2.1 and ^2.2). 
However, there is no such scatter for the mass luminosity 
relation derived from galaxy-galaxy lensing observations, 
as this relation is the result of a global fit to a halo model 
where each galaxy provide little information. Therefore, for 
this relation we assume the same scatter as for the TF 
relation. 

In Figure[9]we show the uncertainty of the magnification 
as a function of the magnification obtained with the TF/FJ 



8 The maximum value of 125" was limited by the maximum 
number of lenses that Q-LET could handle, this has however 
essentially no impact on this analysis: the brightest galaxy of 
our sample placed at 125" of a supernova line of sight causes a 
magnification fi m < 2 10~ 5 . 
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relations. Wc obtain a relative uncertainty on the magnifi- 
cation of 17%. The most important source of uncertainty 
comes from the scatter in the mass-luminosity relation. The 
uncertainties due to redshift are small and represent a rel- 
ative uncertainty of about 5%. 



--.0.04 



O) 

o 




d"0.02 



-0.1 
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Fig. 9. Magnification uncertainty as a function of magnifi- 
cation both expressed in magnitudes. The line is a straight 
line fit to the data. The relative uncertainty is about 17%. 



5. Signature of the lensing effect 

5.1. Signal discriminant 

As a criterion for a lensing signal detection we have chosen 
to calculate the weighted correlation coefficient. 



cov(p mi r) 

III 

)var(r) 



(G) 



with p m = — 2.51og 10 /x, /i being the gravitational magnifi- 
cation factor, and r is the residual to the Hubble diagram 
(hereafter Hubble residual). The weighted covariance of two 
variables x and y, can be written as : 



cov(x, y) 



xy 



(7) 



where w is the weight assigned to each datapoint and x and 
y are the weighted means. Using random lines of sight (ie. 
random source positions in the galaxy catalog), we have 
found that weighting with the inverse of the Hubble resid- 
ual variance is optimal for the signal detection. Note that 
any scaling applied to p, m and hence any common offset ap- 
plied to the mass-luminosity relations will not change the 
weighted correlation coefficient and thus not change the de- 
tection significance. 

5.2. Expected Signal 

Before presenting the results, it is useful to have an idea 
of what to expect. We now have all the tools needed to 
perform detailed Monte Carlo simulations using the true 
galaxy catalog and give precise predictions both for the 
SNLS 3-year data set and the full SNLS sample. 



Magnification distributions for supernovae at different 
redshifts were generated by calculating the magnification 
factor for a sample of simulated supernovae randomly posi- 
tioned in the true galaxy catalog. In this way it was possi- 
ble to simulate the magnification distribution for a sample 
of supernovae with a given redshift distribution, chosen as 
the actual sample redshift distribution. To generate a corre- 
lated sample, i.e assuming a correlation between the Hubble 
residual and the magnification, we first computed a "true" 
magnification using the TF/FJ relations applied to galaxies 
along the line of sight. We then drew an expected magnifi- 
cation using the uncertainty described in 4.4.2 The corre- 
sponding Hubble residuals were eventually generated from 
the expected magnification and a random offset of r.m.s 
0.16 mag in accordance with the SNLS data for which the 
Hubble residuals exhibit an r.m.s of 0.16 mag. An uncor- 
rected sample can be generated by randomly assigning a 
Hubble residual to each expected magnification. 

We first simulated a large number of correlated samples 
similar to the real one with 171 supernovae. We compared 
the distribution of the weighted correlation coefficient of 
these samples with that of uncorrelated samples 



(see Figure 10). We find that for the current sample 
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200 




-0.2 0.2 0.4 

weighted correlation coefficient (p) 



0.6 



Fig. 10. Simulated distribution (using real lines of sight) of 
the weighted correlation coefficients for correlated samples 
(dotted) and uncorrelated samples (solid), for the 3-year 
SNLS sample. There is 50 % probability of finding a 2.5 
a significance correlation or better (p — 0.22) and 35 % 
probability of detecting a 3 a signal (p = 0.25). 



there is 50 % probability of finding a 2.5 a significance 
correlation or better and 35 % probability of detecting a 
3 a signal. 

For the final SNLS sample we expect ~400 spectroscop- 
ically confirmed Type la supernovae and ^200 photometri- 
cally identified Type la supernovae. Due to masking, about 
a third of the supernovae will be rejected. If we do sim- 
ulations for 400 supernovae with the same parameters as 
the current sample we find that there is 80% probability of 
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detecting a 3 a signal or more (see Figure 

Another question to be addressed is how the errors influ- 
ence the possibility of a signal detection. The errors on the 
magnification are already small. Monte Carlo simulations 
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for the simulated 



Fig. 11. Same distributions as Figure 
full SNLS sample. There is 80% probability of detecting a 
3 (7 signal (p = 0.22) or better. 



performed with various scaled errors of the magnification 
support that the scatter in the magnification has 

little impact on the signal detection. The signal detec- 
tion is highly dominated by the scatter in the SN Hubble 
residuals which is not likely to decrease significantly in the 
near future. As a consequence improving the signal detec- 
tion will require better statistics and, if possible, higher 
rcdshift SNe. 
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5.3. Results 

5.3.1. The SNLS supernovae 

The 3-year data release contains 233 spectroscopically con- 
firmed Type la supernovae in the redshift range 0.2-1.05 
after quality cuts. 

The Hubble residual of the supernova is the difference in 
distance modulus of the supernova and the best cosmology 
fit. The distance modulus is estimated from the supernova 
brightness and other parameters resulting from a fit of the 



SALT2 model (Guy et al. 20071 to the supernova light 



curves. The uncertainty of a distance modulus combines 
three sources: the photometric me asurement uncerta inties, 
the light curve model scatter (see Guy et al. (2007) for a 
detailed definition), and the so-called intrinsic scatter which 
expresses our lack of complete physical understanding of 
SN la. This intrinsic scatter is chosen so that the minimum 
X 2 of cosmological fits to the Hubble diagram matches the 
number of degrees of freedom. This ensures that the quoted 
uncertainties of Hubble residuals properly describe their 
actual scatter. 



5.3.2. Magnification 



Figure [12] and 13 show the magnification of each SN as 
a function of redshift and the magnification distribution 
respectively. As expected, most SNe are demagnified with 
respect to a homogeneous universe and some are signifi- 
cantly magnified. Moreover, the magnification distribution 
peaks at a value slightly lower than one and presents a 
long magnification tail. Images of galaxies along the line- 
of-sight of 3 of the most magnified supernovae (04Dliv, 



Fig. 12. Magnification factor of the SNLS supernovae ver- 
sus redshift. Results based on the K06 luminosity-mass re- 
lation (top) and on TF and FJ relations (bottom). Most of 
the SNe are slightly demagnified whereas some are signifi- 
cantly magnified. 



03D4cx and 04D4bq) in the 3-year data set are shown in 
Figure [14] (see also Figure [| & [3} . All 3 SNe have one or 
several massive galaxies very close to the line-of-sight lo- 
cated roughly halfway between the SN and us, causing the 
high magnification. The 8 most magnified supernovae along 
with information on the most important galaxies causing 
the magnification are listed in table [T] 

5.3.3. Correlation 

We are searching for a correlation between the Hubble 
residuals from a best fit cosmological model and the es- 
timated magnifications of the supernovae based on fore- 
ground galaxy modeling. In Figure 15 we show a plot of 
the Hubble residuals of the 171 SNe versus the estimated 
magnification. 

The weighted correlation coefficient for this sample is 
p = 0.12 using the K06 relation and p = 0.18 using the TF 
and FJ relations respectively. To evaluate the strength of 
the correlation we calculate the distribution of the weighted 
correlation coefficient for an uncorrelated sample and com- 
pare it with the obtained value for our sample (see Figure 
16). The uncorrelated samples are drawn by randomly as- 
sociating Hubble residuals and expected magnifications of 
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Table 1. The most magnified supernovae and the characteristics of the galaxies dominating the magnification. 
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It is tempting to attribute the stronger detection of the 
TF/FJ relation to the fact that this method considers sep- 
arately elliptical and spiral galaxies. In order to test the 
influence of the separation, we run the analysis again, but 
with a random galaxy type assignment (however preserving 
the measured proportions). We then find that the signifi- 
cance of the detection drops from 2.3 to 1.4 a, and conclude 
that the higher significance of our detection using TF/FJ 
relations can be attributed to the distinction between spiral 
and elliptical galaxies, rather than chance. So, our primary 
result is the detection of supernovae lensing using TF/FJ 
relations at the 99% CL. Note tha t sim ulations assuming a 
perfect galaxy typing (see section |5.2[) give rise to a mean 
detection level of 2.5 a (see Figure |l0|). The 2.3 a detec- 
tion we find is thus close to the optimal so that further 
improvements on galaxy typing should not give rise to a 
much larger significance of the signal. 

In order to test the scale of galaxy mass estimates, we fit 
the slope a relating Hubble residuals and expected magnifi- 
cations : < r >= a < /i m > and find a = 0.65 ± 0.30. Hence 
the data are consistent with the TF/FJ mass-luminosity re- 
lations at the 1.2 a level, with a precision of 30%. 

Using random lines of sight in the real data, we can 
estimate the increase of Hubble diagram scatter expected 
from gravitational lensing, as a function of redshift. This 
is shown in Figure 17 for the TF/FJ relations, which can 
be roughly describea as cr(^ m ) = 0.08 x z. Alternatively, 
if we use the value of a derived from a fit of the relation 
between Hubble residuals and magnifications, we obtain a 
lower value a(p m ) — (0.05 ± 0.022) x z. 



Fig. 13. Magnification distribution of the SNLS super- 
novae. Results based on the K06 luminosity-mass relation 
(top) and on TF and FJ relations (bottom). The magnifi- 
cation distributions peaks at a value slightly below 1 and 
presents a high magnification tail. 



the real sample. The probability of finding a larger weighted 
correlation coefficient than the measured one from an un- 
corrected sample is 5% using the K06 relation and 1% 
using the TF and FJ relations, corresponding to 1.6 and 
2.3 a detections respectively. 



6. Conclusion 

We have calculated the expected magnification of the SNLS 
3-year sample from the foreground galaxy properties and 
searched for a correlation with residuals from the Hubble 
diagram. A correlation is detected at the 99% CL, compat- 
ible with a slope of 1. 

The expected magnifications cause an extra scatter in 
the Hubble diagram approximated by 0.08 x z from the 
TF/FJ relations we used, which becomes (0.05 ±0.022) x z 
once these TF/FJ relations are calibrated with the super- 
nova data. We show that separating the galaxy sample into 
a blue and a red population based on a U—V color cut in- 
creases the significance of the detection. This is due to the 
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fact that the mass estimate and hence the induced magni- 
fication is significantly different for a red and a blue galaxy 
of same luminosity (the TF and FJ relations). Simulations 
also point to the fact that a signal detection is dominated 
by the number of SNe, their redshift distribution and the 
scatter in the SN residuals. Reducing the scatter in the esti- 
mated magnification by increasing the photometric redshift 




m * * 
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Fig. 14. From the top to the bottom, SN 04Dliv at a red- 
shift of 0.998, SN 03D4cx at a redshift of 0.949 and SN 
04D4bq at a redshift of 0.550. In red, the SN and its host. 
In blue the foreground galaxies labeled according to the 
closest elliptical distance to the SN and in green the back- 
ground galaxies. All images span 30" (horizontally) and are 
oriented north-up/east-left. 
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Fig. 15. The Hubble residuals as a function of magnifica- 
tion for the SNLS. Results based on the K06 luminosity- 
mass relations (top) and on TF and F J relations (bottom) . 
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Fig. 16. Distribution of the weighted correlation coeffi- 
cient for an uncorrelated sample when using TF and FJ re- 
lations (in black). The red line indicates the value p = 0.18 
obtained for our sample. The measured correlation signifi- 
cance is 2.3 a (99% CL). 



precision or reducing the scatter in the TF and FJ relations 
have little effect on the probability of a signal detection. 

Finally, simulations using the true galaxy catalog 
show that using the full SNLS data set (^400 expected 
spectroscopically confirmed Type la SNe and ~200 photo- 
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Fig. 17. Expected magnification scatter from random lines 
of sight as a function of the source redshift, using the TF/FJ 
mass-luminosity relations. 



metrically identified Type la SNe) there is 80% chance of 
detecting a 3 a signal or more. 
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